Low temperature value of the upper critical field H c2 (0) of isotropic single band metals 
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The strong coupling Schossmann-Schachinger theory for the upper critical field is considered me- 
thodically from the user point of view. We check the accuracy of approximate formulas for the upper 
critica l field. In particula r, we explain in detail as a recently proposed convenient expression (Shulga 
et al. 



ond- mat/0103154) can be justified. The connection of H C 2(0) with the critical temperature 



T c is considered and the similarity with the Allen-Dynes T c -formula is shown explicitly. 
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I. INTRODUCTION 

The standard isotropic single band (ISB) Eliashberg 
model jj] is the most developed part of the modern the- 
ory of superconductivity. It describes quantitatively the 
renormalization of the physical properties of metals due 
to the electron-phonon interaction. The input mate- 
rial parameters are the density of states at Ep, N(0), 
the bare Fermi velocity vf, the impurity scattering rate 
"Jimp, the paramagnetic impurity scattering rate y m , the 
Coulomb pseudopotential 0.1), and the electron- 

phonon spectral (Eliashberg) function a 2 F(U). Thereby 
the validity of the Born approximation is assumed, i.e. 
the system is far from the onset of localization. 

These days the material parameters N(0), the aver- 
aged Fermi velocity, and the spectral function become 
available from microscopic calculations just very soon af- 
ter the discovery of a new superconducting compound. 
The standard procedures of analysis of many physical 
properties of metals include these results as input pa- 
rameters. For example, the density of states N(0) is 
important for the description of the specific heat data. 
Similarly, the calculated plasma frequency is used for 
the interpretation of the optical and transport proper- 
ties. As a rule, the forthcoming optical measurements do 
confirm its value later on. At the same time, the upper 
critical field H c2 (0), the fundamental quantity of a type- 
II superconductors is usually analyzed by experimental- 
ists in the framework of the phenomenological Ginzburg- 
Landau (GL) theory which states H c2 (0) = <J>o/27r£ 2 . 
Here $o is the flux quantum and £ is the GL-coherence 
length, the material parameter of the GL theory. 

Within the Eliashberg ISB model H C 2(T) can be com- 
puted using the Schossmann-Schachinger self-consistent 
theory This theory treats H c2 (T) as the field depen- 
dence of the critical temperature T C {H). At H = the 
linear lf C 2(T)-problem is reduced to the equations for T c . 
From a numerical point of view the eigenvalue problems 
for T c and H C 2(T) are identical and can be computed 
by the same code. In practice, however, the linearized 
equations for T c are frequently not exploited even by 
teams who perform tedious numerical band structure or 



Eliashberg-theory calculations. Instead the well-known 
phenomenological Allen-Dynes expression for T c is widely 
used. In our opinion a similar approach could be applied 
also to H C 2(0) calculating it from a simple formula. 

Here, for the sake of clarity and simplicity the Pauli 
limiting, unimportant in many cases of practical interest, 
is ignored. In the present paper we examine five simple 
expressions for H C 2(0) by comparison with the numerical 
solutions of Eqs. pj-pl One of these formulas valid in the 
clean limit is presented in the present work for the first 
time. It is shown that at least three formulas fit reason- 
ably well the numerical values of H C 2 (0) in a wide range 
of coupling strengths A. The deviations do not exceed 
12%. Extensions to cases of arbitrary impurity content 
are considered as well. An application example of the 
approximate phenomenological formula is provided. 



II. GENERAL EQUATIONS AND 
APPROXIMATE SEMI-ANALYTICAL 
FORMULAS 

To take into account fully retardation and strong cou- 
pling effects in calculating H C 2, we have solved numeri- 
cally the linearized equations proposed in Ref. H 
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Here ui n =irT(2n + 1) are the Matsubara frequencies, e 
and A(n) denote the electron charge and the supercon- 
ducting order parameter, fi and c are the Plank's con- 
stant and the speed of light. 
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From the numerical point of view the equations of 
the Werthamer-Helfand-Hohenberg (WHH) || and the 
Schossmann-Schachinger theory have the same order of 
complexity. At the same time the WHH equations al- 
low full analytical solutions, both in the clean and dirty 
limits. In particular, the clean limit WHH orbital upper 
critical field for the ISB model reads 
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where ks and 7=0.577 denote the Boltzmann's and the 
Euler's constants. In practical units this formula can be 
rewritten as 



#™(0)[Tesla] = 



0.0231T 2 [K] 0.02T C 2 [K] 



,[10 7 cm/s] 



,[10 7 cm/s] 
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Within the BCS theory, where 2A = 3.53T C , the formula 
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is identical to the basic Eq. (^|). It was used by Sung j|] 
for an estimation of vp in Nb. 

McMillan and Werthamer || were the first who intro- 
duced the first order strong coupling corrections (Z(0) = 
1 + A) 
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They conclude, that "this produces a substantial over- 
all increase in H c2 from that which would be obtained 
using the bare band Fermi velocity, but being relatively 
temperature-independent the shift cancels out of h(t)." 
These days Eq. (||) is often regarded as a BCS formula 
since the "bare" Fermi velocity has to be replaced by 
the "dressed" Fermi velocity v* F — v F /(l + A) in Eq. Jq. 
Similarly, the renormalized form of Sung's formula rtq) 
reads 
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For the sake of completeness, we note that in the presence 
of magnetic pair-breaking the considered renormalization 
of vf is different 

Since there are no analytical solution of Eqs. ap- 
proximate expressions should be used to estimate H c2 (0) . 
The parameters in such formulas are determined by the 
fit of numerical data. The first expression JlJ] reads 
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where uj\ n is the characteristic phonon frequency intro- 
duced by Kirzhnits et al. and used later on in Ref. 



|T3[ by Allen and Dynes. A similar H c2 (0) expression was 
proposed by Masharov || who obtained however different 
values for the numerical coefficients. 

Our factorable formula, first presented here, 
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is suitable for applications together with other factorable 
expressions. In practical units it reads 



ff,g(0)[Tesla] 
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Note that the choice of the coefficients in Eq. ( |l3|) is not 
unique. For example, the following expression 
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demonstrates even better agreement with the numerical 
data. Eq. [l4] is especially designed for the Q-check 11 in 
the spirit of the Langmann's BCS Eq. 24 in Ref. M. 



III. RESULTS AND DISCUSSION 

To check the accuracy of Eqs. (|9|-^4]), we solved the 
Eq. (|l|j^) employing an Einstein spectral function with 
various coupling constants A = 0.7 — 2.9 and standard 
values for the Coulomb pseudopotential /i*. Then the 
H C 2 (O)-values from the approximate formulas (p[jl4]) were 
normalized by the "exact" numerical values. The result 
is shown in Fig. ^. If the symbol occurs below the unity 
line, the corresponding approximate expression underes- 
timates H C 2(0) and vice versa. The phenomenological 
formulas Eqs. (|Tl[-^4[) fit well the numerical data. The 
deviations from unity do not exceed 12% in a wide range 
of coupling strengths. The formulas Eqs. (|9||lC|) being 
identical within BCS theory deviates differently and ex- 
hibit relatively lower accuracy. 

A detailed discussion and collections of the BCS up- 
per critical field formulas are given for instance in Ref. 
11. Here we discuss only the structure of weak cou- 
pling formulas Eqs^ (|g|-|To|) with respect to the general 
Bergman-Rainer ansatz 



H c2 (0)=H^ i '(0)f(a 2 F(n)^*,T c ). 



(15) 



For the reason of dimensionality H^ 2 CS (0) is proportional 
to the square of some typical energy. It could be T c , the 
gap Ao , the Debye energy ojd , or the average phonon fre- 
quency uji n , and there is a total equivalence between them 
within the BCS theory due to its universal character. 
Within the Eliashberg theory the situation is different. 
Here Eq. (^|) is more preferable than the renormalized 
Sung's formula ( [To| ) (see Fig. |Tj). In fact, at given A the 
jjMcM-whh ^ p 0m t s accumulate at nearly the same 
position for different values of /i* in comparison with the 
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FIG. 1. The upper critical field, normalized on the corre- 
sponding "exact" Schossmann-Schachinger numerical values, 
vs. coupling constant A. Approximate basic WHH formula 
(V, Eq. |9|) ,renormalized Sung formula (A, Eq. ho|), Carbotte 
formula (O, Eq. |ll]), and our formulas (O; Eq. [w[and □, Eq. 
|l3| ) were indicated. The fitted functions (1 + A) -0,231 and 
0.02(1 + A)"°' 4 /0.0231 were shown by solid and dash-dotted 
lines. Identical symbols for the same A correspond to different 
standard values of the Coulomb pseudopotential mu* =0.06, 
0.1, and 0.14, respectively. 



H^ cMS (0) data which spread out. This accumulation is 
a clear advantage of the basic Eq. (^|). A similar choice in 
favor of Eq. (||) was made by Carbotte jl| without com- 
ments. We believe that the weak sensitivity of the ratio 
H C 2(0)/T^ to details of a 2 F(uj) and the actual value of 
/i* is due to the fact that H C 2(0) and T c are the solution 
of the same linearized system of Eqs. ([l|||) jf2j. 

In Sec. U various H C 2(0) formulas are presented in 
different styles. The first group consists of "textbook- 
like" expressions (|o|p| JTTlJT^ ) which contain explicitly fun- 
damental physical and mathematical constants. Eqs. 
0,|l3],[l4| written in practical units are preferred by insid- 
ers working in the field. Sung's formulas j^JlQ|) occupy 
an intermediate position. In our opinion the complete 
"textbook-like" presentation is necessary when a formula 
is derived for the first time from a general theory without 
simplifying additional assumptions. Otherwise, all styles 
are equivalent. In the present context a proof means a 
comparison of values given by a formula and a related 
computer code. 

In view of a reasonable accuracy the question might 
arise: "Should we truncate higher digits in the prefactor 
B x = 0.0231 of Eq. @ to B 2 = 0.02 before substitution 
into Eq. [ll] or not?" because the Fermi velocity is not 
constant across the Fermi surface, and it is usually known 
with some uncertainty. The alternative possibilities led 
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to Eqs. (|12j,|13|) and Eq. 

The Carbotte formula ( |ll[ ) also derived from Eqs. 
(|l[]|) through a series of simplifying assumptions, in- 
cludes second order strong coupling corrections in terms 
of the small parameter T c /un n . The leading linear term 
1.44r c /wi n is introduced on a pure phenomenologically 
basis to get the best fit of numerical data (see remark 
below Eq. (7.15) in Ref. (![]). The Carbotte formula con- 
tains the additional parameter T c /oj\ n in comparison with 
"bare" BCS one and is not factorable. 

For the overwhelming majority of ISB superconductors 
the Allen-Dynes formula 
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describes T c within an accuracy of about 5%. In a rigor- 
ous sense, it can be reduced to the BCS expression 
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at A — /i* — > only, when T c — >0. The coupling constant 
A, the Coulomb pseudopotential /i*, and the characteris- 
tic boson energy lu c = w;„/(1.13* 1.2) are the parameters 
which enter Eqs. [h| [jj]. Despite the formal "contradic- 
tion" to the standard BCS Eq. (|l7|), the Allen-Dynes for- 
mula is widely used, especially for qualitative discussions 
where possible slight uncertainties may be ignored. Our 
formulas Eqs. (0-14) were designed in the same manner 
and can be treated as the counter part of the "Allen- 
Dynes formulas" for the clean limit H C 2(0) of isotropic 
single band metals. 

H c2 (0) does not depend on N(0) and in view of Eq. 
(H), H c2 (0) cx 1/vp. As discussed above, the relation 
H C 2(0) oc T 2 captures with high accuracy the depen- 
dences on /i* and on the shape of the spectral function. 
The remaining material parameter A measures the overall 
strength of the spectral function a 2 F{oj). Since A enters 
the basic BCS formula as the mass renormalization (1+ 
A) it is natural at first to make for the correction function 
entering Eq. (Q5h the following ansatz 



/ = a + Ar 



(18) 



The chi-by-eye fl5f fit is shown in Fig. |] for the prefactor 
values 5i=0.0231 and #2=0.02. The round-off before 
the fit leads to A2 «0.4 in comparison with Ai=0.23 
obtained for BCS B2 ~0.0231. The additional calcula- 
tions using various spectral functions and Coulomb pseu- 
dopotentials fj,* yield a deviation of SA «0.08. Thus we 
adopt Ai =0.2 and A2=0. 4. Notice that an attempt to 
renormalize the Fermi velocities by the (1 + \) 1+A / 2 w 
(1 + (1 + A/2)X) and the reciprocal BCS-coupling con- 
stant as 1.04(1 + A)/[A — /i*(l + 0.62A)] are in a rigorous 
sense incorrect, since the meaning of "dressing" and mass 
renormalization correspond to linear relations. 
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FIG. 2. The calculated upper critical field for pi*=0.13 
and A =1 (□), 2 (o), 3 (O) (Eqs. normalized on its ap- 
proximate clean limit value (Eq. n3l) vs normalized impurity 
scattering rate 7i mp /r e (l + A). Solid line indicates BCS ap- 
proximation (Eq. [19). 



Formally speaking, the discussed above approximate 
formulas should be supplemented by a clean limit cri- 
terion. We remind the reader that there is no natural 
criterion for H c 2 in the weak coupling BCS regime. In- 
stead the general case of arbitrary impurity scattering 
rate can be described approximately by a sum of two 
terms It is shown [lj that the second order cor- 

rections with an accuracy of about 6% coincide in the 
clean and extreme dirty limits. Hence, one expects that 
H C 2(0,jim P )/H^2 ean (0), the upper critical field at finite 
limp normalized on its clean limit value, can be approx- 
imated by the BCS formula. We solved Eqs. E[@ for 
Einstein spectra. The results are shown in Fig. For 
comparison the ratio of the BCS clean and dirty limit 
expressions 
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is depicted too. One realizes that this simple BCS-type 
relationship holds in the moderate strongly coupled case, 
too. 

To summarize, for nearly isotropic single band met- 
als the upper critical field at low temperature can be 
estimated with high accuracy by rather simple semi- 
analytical formulas Eqs. (pl]-p4[) and Eq. (|l9|). 
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APPENDIX A: 



AN INSTRUCTIVE EXAMPLE 
FROM MgB 2 



Recently Manske et al. Q suggested that experi- 
mental data for the novel superconductor MgB2 which 
we interpreted in terms of a multi-component gap ap- 
proach M, can be also described within the standard 



4 



isotropic single-band Eliashberg model using Eq. (11). 
They reported 2A /k B T c w 4.1 and H c2 (0)= 14 T in 
seemingly good agreement with experimental data for 
H c 2(0). The other reported quantities are u^„=30 meV, 
T c =30 K. Using their input (material) parameters pre- 
sented in that paper we checked their calculations and 
find that their model gives actually u>i n =53 meV, T C —2Q 
K, 2A /fc s T c =3.8, and H c2 {T = 0) = 1.2 T in sharp 
contrast with their set mentioned above. Thus, at first, 
we confirm our earlier result |9| that the standard ISB 
model can not be applied to MgB 2 . At second, we note 
that the use of formulas written in practical units might 
be helpful to avoid such incorrect estimations. 
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